- Author:
- Soroush <ssaf006@aucklanduni.ac.nz>
- Date:
- 2021-04-08 15:18:19+12:00
- Desc:
- modified svg
- Permanent Source URI:
- https://models.fieldml.org/workspace/43b/rawfile/994dbbb88a761fe3a62bfaf8352da94103650bb2/BG Tutorial Biochemical Systems/OHb.cellml
<?xml version='1.0' encoding='UTF-8'?>
<model name="Oxyhaemoglobin" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#" xmlns:xlink="http://www.w3.org/1999/xlink">
<import xlink:href="Units.cellml">
<units name="per_mol" units_ref="per_mol"/>
<units name="J_per_mol" units_ref="J_per_mol"/>
<units name="J_per_mol_K" units_ref="J_per_mol_K"/>
<units name="mol_per_s" units_ref="mol_per_s"/>
</import>
<component name="state">
<variable initial_value="0.0" name="t" units="second"/>
<!-- Parameters-->
<variable initial_value="8.324" name="R" units="J_per_mol_K"/>
<variable initial_value="300" name="T" units="kelvin"/>
<variable initial_value="1.0" name="K_q_R_Hb" units="per_mol"/>
<variable initial_value="1.0" name="K_q_R_OHb" units="per_mol"/>
<variable initial_value="1.0" name="K_q_R_OOHb" units="per_mol"/>
<variable initial_value="1.0" name="K_q_R_OOOHb" units="per_mol"/>
<variable initial_value="1.0" name="K_q_R_OOOOHb" units="per_mol"/>
<variable initial_value="1.0" name="K_q_T_Hb" units="per_mol"/>
<variable initial_value="1.0" name="K_q_T_OHb" units="per_mol"/>
<variable initial_value="1.0" name="K_q_T_OOHb" units="per_mol"/>
<variable initial_value="1.0" name="K_q_T_OOOHb" units="per_mol"/>
<variable initial_value="1.0" name="K_q_T_OOOOHb" units="per_mol"/>
<variable initial_value="1.0" name="K_q_RBC_R_O2" units="per_mol"/>
<variable initial_value="1.0" name="K_q_RBC_T_O2" units="per_mol"/>
<variable initial_value="1.0" name="K_Re_RT_Hb" units="mol_per_s"/>
<variable initial_value="0.1" name="K_Re_RT_OHb" units="mol_per_s"/>
<variable initial_value="0.001" name="K_Re_RT_OOHb" units="mol_per_s"/>
<variable initial_value="0.0001" name="K_Re_RT_OOOHb" units="mol_per_s"/>
<variable initial_value="0.00001" name="K_Re_RT_OOOOHb" units="mol_per_s"/>
<variable initial_value="4.0" name="K_Re_R_OHb" units="mol_per_s"/>
<variable initial_value="1.5" name="K_Re_R_OOHb" units="mol_per_s"/>
<variable initial_value="0.33" name="K_Re_R_OOOHb" units="mol_per_s"/>
<variable initial_value="0.25" name="K_Re_R_OOOOHb" units="mol_per_s"/>
<variable initial_value="4.0" name="K_Re_T_OHb" units="mol_per_s"/>
<variable initial_value="1.5" name="K_Re_T_OOHb" units="mol_per_s"/>
<variable initial_value="0.33" name="K_Re_T_OOOHb" units="mol_per_s"/>
<variable initial_value="0.25" name="K_Re_T_OOOOHb" units="mol_per_s"/>
<!-- State variables-->
<variable initial_value="1.0" name="q_R_Hb" units="mole"/>
<variable initial_value="0.0" name="q_R_OHb" units="mole"/>
<variable initial_value="0.0" name="q_R_OOHb" units="mole"/>
<variable initial_value="0.0" name="q_R_OOOHb" units="mole"/>
<variable initial_value="0.0" name="q_R_OOOOHb" units="mole"/>
<variable initial_value="0.0" name="q_T_Hb" units="mole"/>
<variable initial_value="0.0" name="q_T_OHb" units="mole"/>
<variable initial_value="0.0" name="q_T_OOHb" units="mole"/>
<variable initial_value="0.0" name="q_T_OOOHb" units="mole"/>
<variable initial_value="0.0" name="q_T_OOOOHb" units="mole"/>
<variable initial_value="1.0" name="q_RBC_R_O2" units="mole"/>
<variable initial_value="0.0" name="q_RBC_T_O2" units="mole"/>
<variable name="u_R_Hb" units="J_per_mol"/>
<variable name="u_R_OHb" units="J_per_mol"/>
<variable name="u_R_OOHb" units="J_per_mol"/>
<variable name="u_R_OOOHb" units="J_per_mol"/>
<variable name="u_R_OOOOHb" units="J_per_mol"/>
<variable name="u_T_Hb" units="J_per_mol"/>
<variable name="u_T_OHb" units="J_per_mol"/>
<variable name="u_T_OOHb" units="J_per_mol"/>
<variable name="u_T_OOOHb" units="J_per_mol"/>
<variable name="u_T_OOOOHb" units="J_per_mol"/>
<variable name="u_RBC_R_O2" units="J_per_mol"/>
<variable name="u_RBC_T_O2" units="J_per_mol"/>
<variable name="u_R1_OHb" units="J_per_mol"/>
<variable name="u_R1_OOHb" units="J_per_mol"/>
<variable name="u_R1_OOOHb" units="J_per_mol"/>
<variable name="u_R1_OOOOHb" units="J_per_mol"/>
<variable name="u_T1_OHb" units="J_per_mol"/>
<variable name="u_T1_OOHb" units="J_per_mol"/>
<variable name="u_T1_OOOHb" units="J_per_mol"/>
<variable name="u_T1_OOOOHb" units="J_per_mol"/>
<variable name="v_R_OHb" units="mol_per_s"/>
<variable name="v_R_OOHb" units="mol_per_s"/>
<variable name="v_R_OOOHb" units="mol_per_s"/>
<variable name="v_R_OOOOHb" units="mol_per_s"/>
<variable name="v_T_OHb" units="mol_per_s"/>
<variable name="v_T_OOHb" units="mol_per_s"/>
<variable name="v_T_OOOHb" units="mol_per_s"/>
<variable name="v_T_OOOOHb" units="mol_per_s"/>
<variable name="v_RT_Hb" units="mol_per_s"/>
<variable name="v_RT_OHb" units="mol_per_s"/>
<variable name="v_RT_OOHb" units="mol_per_s"/>
<variable name="v_RT_OOOHb" units="mol_per_s"/>
<variable name="v_RT_OOOOHb" units="mol_per_s"/>
<!-- Conservation laws-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>u_R1_OHb</ci>
<apply>
<plus/>
<ci>u_R_Hb</ci>
<ci>u_RBC_R_O2</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>u_R1_OOHb</ci>
<apply>
<plus/>
<ci>u_R_OHb</ci>
<ci>u_RBC_R_O2</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>u_R1_OOOHb</ci>
<apply>
<plus/>
<ci>u_R_OOHb</ci>
<ci>u_RBC_R_O2</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>u_R1_OOOOHb</ci>
<apply>
<plus/>
<ci>u_R_OOOHb</ci>
<ci>u_RBC_R_O2</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>u_T1_OHb</ci>
<apply>
<plus/>
<ci>u_T_Hb</ci>
<ci>u_RBC_T_O2</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>u_T1_OOHb</ci>
<apply>
<plus/>
<ci>u_T_OHb</ci>
<ci>u_RBC_T_O2</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>u_T1_OOOHb</ci>
<apply>
<plus/>
<ci>u_T_OOHb</ci>
<ci>u_RBC_T_O2</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>u_T1_OOOOHb</ci>
<apply>
<plus/>
<ci>u_T_OOOHb</ci>
<ci>u_RBC_T_O2</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_R_Hb</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_R_OHb</ci>
</apply>
<ci>v_RT_Hb</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_R_OHb</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_R_OHb</ci>
<ci>v_R_OOHb</ci>
</apply>
<ci>v_RT_OHb</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_R_OOHb</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_R_OOHb</ci>
<ci>v_R_OOOHb</ci>
</apply>
<ci>v_RT_OOHb</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_R_OOOHb</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_R_OOOHb</ci>
<ci>v_R_OOOOHb</ci>
</apply>
<ci>v_RT_OOOHb</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_R_OOOOHb</ci>
</apply>
<apply>
<minus/>
<ci>v_R_OOOHb</ci>
<ci>v_RT_OOOOHb</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_T_Hb</ci>
</apply>
<apply>
<plus/>
<ci>v_T_OHb</ci>
<ci>v_RT_Hb</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_T_OHb</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_T_OHb</ci>
</apply>
<ci>v_T_OOHb</ci>
<ci>v_RT_OHb</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_T_OOHb</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_T_OOHb</ci>
</apply>
<ci>v_T_OOOHb</ci>
<ci>v_RT_OOHb</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_T_OOOHb</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_T_OOOHb</ci>
</apply>
<ci>v_T_OOOOHb</ci>
<ci>v_RT_OOOHb</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_T_OOOOHb</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_T_OOOHb</ci>
</apply>
<ci>v_RT_OOOOHb</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_RBC_R_O2</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_R_OHb</ci>
</apply>
<ci>v_R_OOHb</ci>
</apply>
<ci>v_R_OOOHb</ci>
</apply>
<ci>v_R_OOOOHb</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_RBC_T_O2</ci>
</apply>
<apply>
<plus/>
<ci>v_T_OHb</ci>
<ci>v_T_OOHb</ci>
<ci>v_T_OOOHb</ci>
<ci>v_T_OOOOHb</ci>
</apply>
</apply>
<!-- Constitutive relations-->
<apply>
<eq/>
<ci>u_R_Hb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_q_R_Hb</ci>
<ci>q_R_Hb</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>u_R_OHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_q_R_OHb</ci>
<ci>q_R_OHb</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>u_R_OOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_q_R_OOHb</ci>
<ci>q_R_OOHb</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>u_R_OOOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_q_R_OOOHb</ci>
<ci>q_R_OOOHb</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>u_R_OOOOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_q_R_OOOOHb</ci>
<ci>q_R_OOOOHb</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>u_T_Hb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_q_T_Hb</ci>
<ci>q_T_Hb</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>u_T_OHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_q_T_OHb</ci>
<ci>q_T_OHb</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>u_T_OOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_q_T_OOHb</ci>
<ci>q_T_OOHb</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>u_T_OOOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_q_T_OOOHb</ci>
<ci>q_T_OOOHb</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>u_T_OOOOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_q_T_OOOOHb</ci>
<ci>q_T_OOOOHb</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>u_RBC_R_O2</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_q_RBC_R_O2</ci>
<ci>q_RBC_R_O2</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>u_RBC_T_O2</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_q_RBC_T_O2</ci>
<ci>q_RBC_T_O2</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_R_OHb</ci>
<apply>
<times/>
<ci>K_Re_R_OHb</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_R1_OHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_R_OHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_R_OOHb</ci>
<apply>
<times/>
<ci>K_Re_R_OOHb</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_R1_OOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_R_OOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_R_OOOHb</ci>
<apply>
<times/>
<ci>K_Re_R_OOOHb</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_R1_OOOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_R_OOOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_R_OOOOHb</ci>
<apply>
<times/>
<ci>K_Re_R_OOOOHb</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_R1_OOOOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_R_OOOOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_T_OHb</ci>
<apply>
<times/>
<ci>K_Re_T_OHb</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_T_OHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_T1_OHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_T_OOHb</ci>
<apply>
<times/>
<ci>K_Re_T_OOHb</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_T_OOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_T1_OOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_T_OOOHb</ci>
<apply>
<times/>
<ci>K_Re_T_OOOHb</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_T_OOOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_T1_OOOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_T_OOOOHb</ci>
<apply>
<times/>
<ci>K_Re_T_OOOOHb</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_T_OOOOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_T1_OOOOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_RT_Hb</ci>
<apply>
<times/>
<ci>K_Re_RT_Hb</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_R_Hb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_T_Hb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_RT_OHb</ci>
<apply>
<times/>
<ci>K_Re_RT_OHb</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_R_OHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_T_OHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_RT_OOHb</ci>
<apply>
<times/>
<ci>K_Re_RT_OOHb</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_R_OOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_T_OOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_RT_OOOHb</ci>
<apply>
<times/>
<ci>K_Re_RT_OOOHb</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_R_OOOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_T_OOOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_RT_OOOOHb</ci>
<apply>
<times/>
<ci>K_Re_RT_OOOOHb</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_R_OOOOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>u_T_OOOOHb</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
</model>